查看原文
其他

手把手教你EEG脑电数据预处理-操作篇

亚内士多德 脑机接口社区 2022-04-26

EEG脑电数据预处理-操作篇

关于EEG数据预处理的原理,前面已经介绍过了,大家可以查看《手把手教你EEG脑电数据预处理-原理篇》。下面是亚内士多德分享的操作篇。


基本步骤

  • 定位通道位置

  • 删除无用通道

  • 滤波

  • 分段

  • 基线校正

  • 重参考

  • 降低采样率

  • 插值坏导

  • 独立主成分分析

  • 剔除坏段



导入数据

选择要导入的数据格式


接下来会弹出一个对话框,是否要对数据进行选择性导入,一般都是全部导入,直接点OK即可。


接下来又会再弹出一个对话框,是否要对数据进行命名。基本上在EEGLAB的每一步操作之后都会弹出这样的对话框,询问你是否需要对新产生的数据进行命名,根据自己需要选择即可。


对数据进行初步认识


  • channel per frame:64   指导入的数据有64个通道

  • frames per epoch:一段数据的总长度,是439880采样点

  • epochs:指当前数据的段数,原始数据还没有进行分段,因此显示只有一段

  • events:检测到当前数据一共有202个events

  • sampling rate(HZ):数据的采样率为1000HZ

  • epoch start和epoch end: 这个的分段是从0秒开始,到439.879秒结束。还没有进行分段所以看这个数值没有意义

  • reference: 指数据的参考点,重参考后会显示重参考的电极点,或者average,目前还没有进行重参考所以是unknown

  • channel locations:是否有对通道进行定位,目前显示没有,定位后会显示为yes

  • ICA weights:是否对数据进行了ICA独立主成分分析,分析后会显示yes

  • dataset size:数据的大小


电极定位

默认文件是standard-10-5-cap385.elp,点击确定即可

channel locations


channel locations


点击OK之后可以进入下一步查看数据

channel locations info


如果检查到某个电极点的坐标信息为空白(如上图所示),可能是该通道的名称在加载的模板文件中找不到,那么

  • 如果确定该通道是根据国际排布系统来确定的,那可以修改label名称后再点击下方的look up locs来更新通道位置信息,比如上图的VEO改为VEOG,然后更新位置信息之后就可以看到了

  • 如果不是国际排布系统,可以向厂家要电极的坐标位置信息

  • 如果知道该电极点的精确坐标,可以手动输入

  • 双侧乳突的电极点位置可能会根据脑电系统的不同而不同,比如TP9 TP10,A1 A2, M1 M2这三对都是指双侧乳突

更新后的电极点坐标


移除无用电极

此处我们要移除掉两个眼电数据

select data


select data


select data


重参考

此处要将数据转换为双侧乳突平均参考,即TP9,TP10

re-reference


re-reference


重采样

不一定每个数据都需要进行重采样,可根据需要进行。

re-sampling


re-sampling


滤波

此处选择默认的滤波器进行0.1-30Hz的滤波

filter


依次在第一框中填入0.1和第二个框中填入30进行计算,两个一起填代表进行带通滤波,建议分开进行

filter


分段和基线校正

extract epochs


extract epochs


点击OK后会弹出下一个要求基线校正的对话框,默认是event前的数据作为基线来进行校正,点击OK即可。

baseline correction


插值坏导

对数据进行检查,如果发现某个通道的数据坏了,可以用插值的方式来进行校正,此处要校正的通道为FC3,选择插值的通道为它周围的四个电极点,F3,FC1,C3,FC5


  • 第一种方法是代码操作

EEG.data(a,:,:) = mean(EEG.data([b c d],:,:)) ;

其中a代表坏掉的电极点数值,b c d 代表要进行插值平均的电极点数值
如果不确定通道数值是多少,可以通过Plot - channel locations - by name 来进行查看。


画出通道图之后,可以通过单击通道名字来查看通道数值

channel name


channel name


所以这句代码为

EEG.data(17,:,:) = mean(EEG.data([8 18 26 16],:,:)) ;

这句代码的意思是指,第17个通道的数值,用第8,18,26,16这四个通道数值的平均值来代替。


  • 第二种方法是面板操作,用eeglab的默认算法进行

interpolate electrodes


interpolate electrodes


ICA

使用ICA算法剔除伪迹,此步操作耗时较长

run ICA


因为我们有一个坏通道的数据是用其他通道的数据生成的,相当于我们损失了一个有效通道数据,而ICA默认是在所有通道都是有效数据的情况下进行的,所以此处我们需要进行一些修改,在'extended', 1后面加上'pca', 59,中间用逗号和空格隔开,逗号要用英文输入法下的

run ICA


点击确定之后会出现一个对话框,显示正在进行ICA计算,计算结束后该对话框会消失

run ICA


在跑完ICA之后建议保存当前数据,在进行伪迹剔除操作。这样做可以保证在剔除不当的时候随时回到未剔除的状态下重新进行选择。

此时会保存为eeglab的.set的格式文件

save dataset


在跑完ICA之后就可以用ICA进行伪迹剔除了

reject ICA


选择要画出的成分数,此处我们选择默认状态

plot ICA


eeglab提示我们一张画布只能画下35个图,剩下的会在第二张画布中画出,点击OK即可

plot ICA


plot ICA


此时我们可以对所有的ICA成分进行查看和标记,点击成分数字,会出现该成分的详细情况。如果觉得该成分代表伪迹成分,想要剔除的话,可以先将它标记起来。标记的方式是单机下方绿色的ACCEPT,单击之后会变成红色的REJECT

plot ICA


这一步需要对每一个成分进行查看和判断,然后将想要剔除的成分先标记起来。ICA可以进行眼电伪迹,肌电伪迹,坏导伪迹以及心电伪迹等的辨识与剔除,但是此处我们只进行眼电伪迹的剔除,我们将第一个成分和第13个成分标记起来,标记了的成分数会变成红色。然后点OK。此时还没有进行剔除操作。

mark ICA


除了上面那种方式之外,还可以使用逐个成分画出的方式来查看和标记

plot ICA


plot ICA



剔除ICA成分


remove ICA


eeglab默认两个框都是一样的信息,此时记得把第二个框的信息给清除掉

remove ICA


在接下来的对话框里,我们可以对剔除了ICA成分和未剔除的数据进行比较

check data


红色代表未剔除的数据,蓝色代表剔除后的数据

所有通道的ERP图


从红色框中的部分可以看到,我们还是很有效地剔除了眼动伪迹的

所有通道单个trial的图


检查确认无误后点击accept进行剔除操作


剔除坏段

  • 绝对阈值法
    此处我们选择剔除掉超出±100μV的trial。

reject epoch


reject epoch


目视检查法


channel data


逐个trial进行查看,如果发现有伪迹过大的trial想要剔除的话,单击该trial,它会变成黄色。选择完所有trials后,点击下方的reject即可。此处选中的两个trials仅仅作为演示用。

reject trials


保存数据

至此EEG数据的预处理就基本完成了,我们可以再次将这个数据保存起来,以供下一步操作。

来源:简书ID:亚内士多德

不用于商业行为,转载请联系后台

若有侵权,请后台留言,管理员即时删侵!

更多阅读

上海交大和MIT提出的软性机械手,可提供实时的触觉控制


手把手教你EEG脑电数据预处理-原理篇

CMU贺斌教授团队提出FAST-IRES技术

学习神经科学,发展人工智能

厦大教授Science发表“仿生纳流离子学”前瞻性论文!

脑电数据的Epoching处理

注意力缺陷障碍可以通过训练大脑来管理

脑机产业迎来“新标准”,CESI发布《脑机接口标准化白皮书》

2021世界机器人大赛脑机接口系列视频讲座之SSVEP相关

投稿通道

如何让你的工作让更多人知晓和受益?

脑机接口社区就是这样一个连接学界、

企业界和爱好者的平台渠道。


区鼓励高校实验室、企业或个人在我们平台上分享优质内容。


稿件要求

稿件系个人原创作品,若已在其他平台发表,请明确标注。

稿件一经录取,便提供稿费!

投稿通道

微信扫码,备注:投稿+姓名+单位

微信交流群,请扫码上方微信

(备注:姓名+单位+专业/领域行业)

QQ交流群:913607986

你的每一次在看,我都很在意!

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存